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Coherence properties of Bose-Einstein condensates offer the potential for improved interferometric 
phase contrast. However, decoherence effects due to the mean- field interaction shorten the coherence 
time, thus limiting potential sensitivity. In this work, we demonstrate increased coherence times 
with number squeezed states in an optical lattice using the decay of Bloch oscillations to probe the 
coherence time. We extend coherence times by a factor of 2 over those expected with coherent state 
EEC interferometry. We observe quantitative agreement with theory both for the degree of initial 
number squeezing as well as for prolonged coherence times. 



o 
o 

(N 
Oh! 

<N 



> 
OS 
O 
O 
OS 
O 

o 

> 

:3 



X 



Experimental requirements for precision atom interfer- 
ometry are well suited to many of the coherence proper- 
ties of Bose-Einstein condensates [H d 0, B 1^- BECs 
possess narrower momentum distributions than those of 
ultra-cold atomic gases, which removes the need for ve- 
locity selection during initial state preparation. The 
longer coherence length of a condensate improves phase 
contrast, and colder temperatures reduce ensemble ex- 
pansion during long interferometer interrogation times. 
Furthermore, for confined atom-interferometers 0, 
which require spatial separation of a wavepacket in close 
proximity to a guiding surface, the superfluid proper- 
ties of a EEC offer an additional advantage. The mean- 
field interaction energy in EECs provides an energy gap 
to external excitations, effectively decoupling the atomic 
proof-mass from the physical sensor. 

On the other hand, the coherence time for EEC in- 
terferometry can be significantly reduced with respect to 
cold atom sources. Prior to separation, two linked con- 
densates have relative number fluctuations which support 
a well defined relative phase. However, when separated, 
the interplay of a large on-site mean-field interaction with 
large number variance causes rapid dephasing 9]. This 
concern has been addressed previously by using either di- 
lute condensates lOj or alternatively, Fermi gases which 
do not suffer from density broadening mechanisms 
It is possible, however, to retain some of the benefits of 
EEC interferometry while minimizing mean-field induced 
decoherence. The generation of atom-number squeezed 
states from a EEC in an optical lattice 12, 13] has offered 
the possibility to create states with reduced sensitivity to 
mean- field decay mechanisms. 

In this work, we study the characteristic time scale 
for which an array of EECs preserves relative phase in- 
formation after becoming fragmented, and we observe 
prolonged coherence times for number squeezed states. 
The coherence time is probed through the decay time of 
a Eloch oscillation, and we find quantitative agreement 
with theoretical predictions. 

The theoretical treatment for a EEC in an optical lat- 
tice begins with the Eose-Hubbard Hamiltonian [l^ . 
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where g(3 represents the mean- field interaction energy 
and 7 is the inter-well tunneling matrix element with 
both terms dependent on lattice depth. Si denotes the 
external potential term and di and di^ represent single 
particle annihilation and creation operators respectively 
at the ith lattice site. Ni denotes local atom number 
which is replaced by the central well occupation N in the 
discussion below. 

In the regime of low lattice potential and correspond- 
ingly large tunneling {Nj ^ gP) the many-body ground 
state of the system is described by a superfluid state 
with local atom number uncertainty cr{N) = \/N . As 
the potential barrier is adiabatically raised, the interplay 
of the interaction and tunneling terms renders number 
fluctuations energetically unfavorable. In the Eogoliubov 
limit, number fluctuations decrease with increasing lat- 
tice potential as crs{N) = ( with com- 
mensurately increasing on-site phase fluctuations. In the 
limit of A^7 <C the system enters the Mott-Insulating 
regime, where the wavefunction at each site approximates 
an atom-number Fock state 13]. 

We fragment this array by frustrating tunneling be- 
tween adjacent lattice sites. In previous work, this has 
been achieved by diabatically raising a large potential 
barrier In this work, we frustrate tunneling by 

applying a large energy gradient across the array. This 
method is preferred since the on-site mean field energy 
(the third term on the right of Eq. 1) is unaffected by 
this process. This sudden localization is analogous to a 
beam splitter which instantly separates an ensemble into 
two distinct paths. The many-body ground state does 
not have time to react to the perturbation, and therefore 
an array of independent, localized many- atom states is 
formed. The initial array state is no longer in the ground 
state, and therefore, the array phase collapses. 

Relative phase dispersion between adjacent wells is in- 
tuitively understood by considering an array initially pre- 
pared in a shallow lattice potential. After fragmentation. 
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each coherent state \a) can be expanded as a superposi- 

tion of atom-number Fock states, \a) = ^ — |n), 
where the phase of each term in the superposition evolves 
as A(j)n^i = ii{n^£i)t/fi^ and ji = ng[3 -\- Si is the lo- 
cal chemical potential. Each number term in the su- 
perposition evolves at a different rate, leading to rela- 
tive phase dispersion with a characteristic decay time Tc. 
The decay of array coherence is seen in the time depen- 
dence of the expectation value of the order parameter 
I (a) 1= ^e-^(^^)'^V2 where = {VNgP)-\9]. This 
coherence time can be prolonged, however, by reducing 
the number of terms in the initial atom-number state su- 
perposition. For squeezed states with reduced number 
uncertainty as{N), the coherence time increases as 
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We measure Tc for different initial number variances 
by studying the decay of coherent Bloch oscillations 
flTl llM [13 . An energy gradient Si = Ei is applied to the 
array where E is written in units of energy. This drives 
an oscillatory response in the quasimomentum, of the 
atomic Bloch state with a period T = h/27TE. Although 
Bloch oscillations are traditionally observed for condi- 
tions where the array is described by bands delocalized 
spatially over many lattice sites, they also occur for spa- 
tially localized wavefunctions described in the Wannier- 
Stark basis 20] . We isolate the lattice sites when E ^ j 
but ensure that E is not too large as to cause particle 
loss through Zener tunneling |^ . The Bloch oscillation is 
observed experimentally by interferometrically following 
the evolution of the relative phase between adjacent wells 
A$ = q\/2?i^ where A/2 is the lattice period. Dispersion 
in momentum space provides a quantitative indication of 
dephasing. 

The apparatus used in this experiment has been de- 
scribed in detail We load 10^ ^^Rb atoms into a 
time-orbiting potential (TOP) trap. Evaporative cooling 
generates a BEG in the |F = 2, rriF = 2) state with 1500 
atoms, density p ~ 10^^ cm~^ and temperature ~ 0.2 
Tc. Atom number is determined with absorptive imag- 
ing with 20% shot-to-shot fluctuations and is consistent 
with the observed condensate fraction as a function of 
temperature. 

We trap the condensate in a 1-D, vertically oriented 
optical lattice. The lattice light is red detuned from the 
^^Rb resonance and has 1/e radii of 60 /im. The poten- 
tial depth, [/, (measured in Er^ where Er = Ti^k'^ /2m 
and k = 27t/X with A = 852 nm) is calibrated using three 
independent methods which all agree to within 10%. We 
first measure U by driving an even parity parametric ex- 
citation from the lowest energy band [2l| . We also mea- 
sure the period of Kapitza-Dirac diffraction by suddenly 
turning on the optical lattice ^22]. Finally, we measure 
the frequency of Josephson tunneling, which is valid for 
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FIG. 1: 

A) Lattice intensity is shown for the experimental sequence. 



The lattice intensity is ramped up m 1 r; 



350 ms and 



then held constant for TbiocH during which time a magnetic 
field gradient is applied. The lattice and magnetic fields are 
turned off and the atoms ballistically expand for a 12 ms Time 
of Flight (TOF) before being imaged with a probe pulse. B) 
Absorptive images indicating a Bloch oscillation are shown. 
C) Peak width vs. energy offset E is shown with widths con- 
verted to units of 21%k. Insets show absorptive images both 
with and without dephasing 



low lattice depths, U < 12 Er 23]. For lattice depths 
explored in this work {5 < U < 24 Er) we calculate 
27r X 4 < 7/n < 27r X 250 Hz, 27r x 0.6 < gp/h < 27r x 1.8 
Hz and 103 < N < 150, and the vertical 1/e condensate 
array radius ranges between 7—10 lattice sites. 

After state preparation in the optical lattice (Fig.QJV), 
we drive a Bloch oscillation by applying a magnetic field 
gradient along the array axis. We probe the coherence 
of the Bloch oscillation by releasing the array, rapidly 
switching off the lattice within 500 ns. The interferomet- 
ric signal is absorptively imaged with single atom detec- 
tion sensitivity. Fig. shows images depicting a Bloch 
oscillation with T = 1.1 ms, taken with U = 10 Er and 
E/n = 27r X 900 Hz. 

We ensure that adjacent sites are decoupled during 
Teioch by varying E and measuring the width of the 
central interference peak. An increased width indicates 
dephasing of the Bloch oscillation. We apply the field 
gradient for 40 ms before releasing the array. We see 
in Fig. that for small E the Bloch oscillation shows 
minimal dephasing. However, for large energy offsets the 
peak width increases, saturating at our resolution limit 
for 7 - where j/h = 27r x 39 Hz for U = 12 Er. 

We next quantify the degree of number squeezing gen- 
erated at a given lattice depth. We adiabatically increase 
the lattice intensity and then interfere the array, generat- 
ing an interference pattern with sharp peaks on top of a 
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FIG. 2: 

A) The fraction of atoms observed in the incoherent back- 
ground of the interference signature is plotted vs. lattice 
depth. This fraction from the simulated interference of an 
array of Gaussian wavepackets is shown with a solid line. The 
calculated quantum depletion vs. lattice depth is shown with 
a dashed line. B) Peak width is plotted vs. Teioch for [/ = 10 
Er and E/h — 27r x 900 Hz. Solid line is a fit to data to 
extract Tc = 14.2 ±1.3 ms. C) Tc is plotted vs. lattice depth. 
Solid (dashed) line denotes theoretical dephasing of number 
squeezed (coherent) states. 



w{t) = Wf — {wf — WQ)e~^^^^''^'^ Wf is the maximum 
observed width representing a fully dephased signature, 
and wq is the peak width prior to any phase dispersion. 

Fig. [2p shows the summary graph plotting Tc vs. 
lattice depth, where increased lattice depth reflects in- 
creased number squeezing. The theoretical Tc for isolated 
coherent states and number squeezed states (Eq. 2) are 
shown with dashed and solid lines respectively. For very 
low lattice depths Tc is longer then expected for either 
coherent state or squeezed state dephasing. This is likely 
due to insufficient isolation between wells when 7 is large. 
However, we cannot significantly increase E at this lattice 
depth without introducing Zener tunneling losses. For 
intermediate lattice depths with near-Poissonian number 
variance, we observe good correlation of Tc with that ex- 
pected for isolated coherent state condensates. 

For deeper lattice depths, however, we measure long 
coherence times which are in quantitative agreement with 
theoretical number squeezed state dephasing. For num- 
ber squeezed states prepared at U = 22.5 Er^ Tc = 
19.3 ± 3.5 ms. This represents an increase of a factor 
of 2.1 over the expected decay time of an array of coher- 
ent states in the same lattice potential. It is interesting 
to note that here squeezing extends the coherence time; 
typically, the enhanced fragility of squeezed states to loss 
mechanisms results in reduced coherence times ^25j . 



broad incoherent background. We extract this incoherent 
fraction by fitting the signal to a function with three nar- 
row Gaussian peaks and a fourth broad peak. For very 
low lattice depths the incoherent background is within 
our noise floor. The observed incoherent fraction shown 
in Fig.EK is a quantitative metric of number squeezing as 
seen by comparison with two theoretical models. First we 
simulate the interference pattern from an array of Gaus- 
sian BEG wavefunctions with ^{N) determined by our 
experimental conditions 12,^. We extract the simu- 
lated incoherent background (shown with a solid line) by 
using the same fit function as used with our experimental 
data. Second, we calculate the degree of quantum deple- 
tion 115], the fraction of atoms not characterized by an 
order parameter, expected at a given lattice depth. We 
equate this depletion fraction with the incoherent frac- 
tion of an interference signal, shown with a dashed line. 
Note that our observed quantum depletion is significantly 
higher than that observed in Ref. 24] due to our densities 
which are nearly two orders of magnitude lower. 

With demonstrated control over number squeezing and 
site localization, we explore the dependence of Tc on ini- 
tial number variance. Using the experimental sequence 
in Fig. we take an absorptive image of the inter- 
fering atoms for different lattice depths with E kept 
constant. We fit the central vertical peak to a Gaus- 
sian to determine the width as a function of TBioch- 
We measure Tc by fitting the data as in Fig. [2^ to 
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FIG. 3: 

A) Peak width vs. u;RF/27r. Tc = 1.6 MHz for a BEG in the 
bare harmonic trap. B) Peak width vs. Tramp 

To eliminate other potential sources of dephasing, we 
investigate the effects of finite temperature and the adi- 
abaticity of our lattice ramp on the interferometric peak 
width observed in Bloch oscillations. We prepare con- 
densates at different temperatures by varying the final 
frequency cjrf of the evaporative cooling stage. As seen 
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in Fig. iniA., we observe no change in peak width for 
27r X 1.40 < cjRF < 27r X 1.50 MHz with U = 16 Er. 
However, at cjrf = 27r x 1.51 MHz we observe a sudden 
onset of phase broadening. To avoid this thermal de- 
phasing regime all data was taken with cjrf = 27r x 1.44 
MHz. Note that this critical temperature in the lattice is 
different from the BEC transition temperature in a bare 
harmonic trap (corresponding to cjrf = 27r x 1.6 MHz). 

In Fig. inp, we investigate the dependence of peak 
width on lattice intensity ramp speed. A balance is re- 
quired to avoid losses due to lattice heating with very 
long ramp times and non-adiabaticity effects with short 
ramp times. We find that peak width is insensitive to 
ramp speed for 80 < Tramp < 620 ms. 
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FIG. 4: 

Peak width is plotted vs. continued hold time in lattice af- 
ter removing the external field gradient, Trephase- Insets are 
absorptive images showing initial dephasing of the phase con- 
trast and subsequent revival. 

As a final consideration, we explore the possibility 
for coherence restoration after complete dephasing. We 
use the same experimental sequence as in Fig. with 
[/ = 10 Er and TBioch = 80 ms, ensuring that dephasing 
has occurred. This time, however, after turning off the 
magnetic field gradient we continue to hold the atoms 
in the optical potential before interfering them. We see 
phase contrast return nearly completely after a rephasing 
time Tr ~ 10 ms as shown in Fig. ^ While the details of 
this rephasing mechanism require further investigation, 
a two-well model predicts this time to be determined by 
the generalized Josephson frequency uoj = y/NgPj. Our 
observed rephasing time is in excellent agreement with 
this prediction for our parameters 26]. 

In conclusion, we have demonstrated that number 
squeezed states in an optical lattice can extend coherence 
times for interferometry. We have obtained quantitative 



agreement with theory both in calibrating the number 
variance of our initial state as well as in observed coher- 
ence times with number squeezed states. For sensitiv- 
ity considerations in practical interferometers, however, 
note that extended coherence times come with the price 
of increased phase noise due to number squeezing. Fu- 
ture work will explore this tradeoff in optimizing absolute 
sensor performance. 
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